Perturbation Theory for Traveling Droplets 
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Motion of chemically driven droplets is analyzed by applying a solvability condition of perturbed 
hydrodynamic equations affected by the adsorbate concentration. Conditions for traveling bifurca- 
tion analogous to a similar transition in activator-inhibitor systems are obtained. It is shown that 
interaction of droplets leads to either scattering of mobile droplets or formation of regular patterns, 
respectively, at low or high adsorbate diffusivity. The same method is applied to droplets running 
on growing terrace edges during surface freezing. 



I. INTRODUCTION 

Spontaneous motion of droplets on solid substrate 
driven by chemical reactions that influence wetting prop- 
erties has been observed in a number of experiments 
[l|, H, S lj| ■ More recently, a reversible setup allowing for 
restoration of substrate properties and thereby making 
permanent motion possible has been realized experimen- 
tally 1^1^. Spontaneous motion has been observed also 
in processes of surface freezing and melting l7i| . The cause 
of chemically driven motion is deposition of a poorly wet- 
table material on the substrate or, alternatively, dissolu- 
tion of an adsorbed wettable layer beneath the droplet. 
As a result, the droplet tends to migrate to the area with 
more favorable wetting properties outside its footprint. 
The direction of motion, originally chosen at random, 
is sustained because the area left behind is either per- 
manently left nonwettable or takes time being restored. 
Chemically driven motion is of interest both in microflu- 
idics applications and as a possible mechanism for the 
formation of surface patterns. 

The theoretical model of chemically driven droplet mo- 
tion presented recently 0, Q combined hydrodynamic 
equations in lubrication approximation with a linear 
reaction-diffusion equation for adsorbed species. The 
model included a precursor layer that served to resolve 
the contact line singularity. The model equations were 
solved numerically, mapping the regimes of motion and 
its dependence on various parameters of the problem. 
The results showed, in particular, that droplets become 
immobile when diffusion and substrate modification are 
sufficiently fast. The computations were, however, re- 
stricted to moderate ratios of the bulk droplet size to the 
precursor layer thickness, since numerical routines are all 
but impossible to implement for realistic very large ra- 
tios. 

In this communication, I present an analytical solu- 
tion of the problem combining an integral solution of 
the reaction-diffusion equation for a steadily propagating 
droplet [I3| with velocity computation using the solvabil- 
ity condition of perturbed hydrodynamic equations. The 
hydrodynamic problem is considered in Section^ The 
approach, based on the lubrication approximation, is the 
same as in the theory of creeping motion of droplets un- 



der the action of externally imposed gradients. The per- 
turbation approach to the problem of motion driven by 
surface inhomogeneities was pioneered by Greenspan [ij 
who had, however, to rely on a phenomenological relation 
for the motion of the contact line. Brochard 12] de- 
rived droplet translation velocity using integral balances 
between various driving forces and viscous dissipation. 
The latter poses a special problem in view of the notori- 
ous contact line singularity, which can be resolved either 
by introducing a slip lengt h 13. .1^ or by allowing for an 
ultrathin precursor layer [13. The latter approach has 
been applied to derive integral conditions for a droplet 
driven either by an external field (gravity) or changes in 
the precursor thickness due to droplet interactions [iq . 
The result, applied later also to 2D droplets driven by 
a difference of advancing and receding contact angles, 
turned out to give even better approximation than more 
elaborate theory using a precise solution of lubrication 
equations to compute the shape of the bulk droplet [13 ■ 

The integral condition for droplet motion will be de- 
rived in Section [11 Bl in a more formal way as a solvability 
condition of perturbation equations, using an eigenfunc- 
tion of the adjoint problem introduced earlier for a 2D 
problem [l^. The contact line singularity will be re- 
solved, however, in a more traditional way through in- 
troducing a slip length, and the solvability condition will 
be re-derived for this model. Both approaches to elimi- 
nating the singularity are indistinguishable on a macro- 
scopic level, leading to a model-dependent logarithmic 
factor 19| , and therefore the result does not need to rest 
on the existence of a macroscopic precursor or depend 
on a precise way the hydrodynamic equations have to be 
modified in the immediate vicinity of the substrate. 

In Sections IIIII and | IV| I shall concentrate on the re- 
versible setup of Refs. 5, 6]. The basic approach is de- 
scribed in Sectional followed by the velocity computa- 
tion for a single steadily propagating droplet both in fast 
(Section IIIIBp and slow (Section IIIIC|I diffusion limits 
and in a general case (Section IIIID|) . The principal re- 
sult obtained in the fast diffusion (slow velocity) limit is 
the existence of a supercritical traveling bifurcation anal- 
ogous to a similar transition in activator-inhibitor sys- 
tems [23, |21| . I further investigate droplet interactions 
on either side of this transition, resulting in relaxation 



to a regular stationary pattern sustained by long-range 
repulsion f Sect ion llV A() or scattering of mobile droplets 
(^Section liyBj) . The problem of motion driven by surface 
freezing or melting is briefly considered in Section Ivl 



be obtained from the first-order equation obtained by ex- 
panding Eq. Q in d. It turns out, however, that a rela- 
tion between the velocity and the contact angle distortion 
can be obtained without actually solving this equation; 
it is sufficient to compute its solvability condition. 



II. SOLUTION METHOD 

A. Lubrication Equations 

The droplet shape is described in the lubrication ap- 
proximation by the thin film equation 



rj 



(1) 



Here gravity and other external forces are neglected, 7 
is the surface tension of the droplet interface, and V is 
the 2D gradient operator in the plane of the substrate. 
The simplest suitable expression for the effective mobility 
function q{h)^ obtained assuming the viscosity 77 of the 
droplet to be much larger than that of the surrounding 
fluid and applying the Navier slip boundary condition, is 



q{h) = ^ ih + 3X), 



(2) 



where A is the slip length. 

The boundary condition on the droplet contour F, i.e. 
the contact line, is n ■ V/i = —9, where n is the outer 
normal to F and 9 is the contact angle. The droplet is 
stationary when the equilibrium contact angle does not 
depend on position explicitly. An asymmetry of the con- 
tact angle caused by the substrate modification sets the 
droplet into motion. We shall assume that the asymme- 
try is weak, so that the motion is slow and the change of 
the droplet shape can be viewed as a small correction. 

Taking the direction of motion as the x axis, we rewrite 
Eq. (^3) in the comoving frame in the dimensionless form 



Sx-Vh^V ■ [q{h)W^h] . 



(3) 



where 5 — Uri/j is the capillary number based on the 
droplet velocity U, as yet unknown, and x is the unit 
vector in the direction of motion. The length scale in 
this equation remains arbitrary. 
Assuming (5 ^ 1, we expand 



h = ho + Shi + ■ 



= -n • V (/lo + (5/ii + . . .)r • (4) 



The zero order function is the stationary solution 
which verifies the Laplace equation V^/iq = 0. The 
solution with a constant contact angle 6*0 is just a 
paraboloidal cap 
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(5) 



where r is the radial coordinate and R is the droplet ra- 
dius. The perturbed shape of the moving droplet should 



B. Translational Solvability Condition 

The first-order equation has a general form 
£/ii + *(a;) = 0, 
which contains the linear operator 

Chi = -V • [g(/io)VV2/ii] 
and the inhomogeneity 



*(a;) 



V/ln 



(6) 



(7) 



(8) 



The operator C is not self-adjoint. The adjoint equation 
defining the translational Goldstone mode (p is 



This equation is verified by the eigenfunction 

ho 



f 



q{ho. 



■ dx. 



(9) 



(10) 



The integration can be carried out along an arbitrary axis 
X, which can be chosen to coincide with the direction of 
motion. Integrating along the two Cartesian axes gives 
two Goldstone modes corresponding to two translational 
degrees of freedom in the plane. 

The solvability condition of Eq. © defining the trans- 
lation speed is obtained by multiplying it by cp and inte- 
grating over the droplet footprint TZ bounded by a con- 
tour F. Since the solvability condition is evaluated in a 
finite region, it includes both the area and contour inte- 
grals. The area integral stemming from the inhomogene- 
ity is evaluated using integration by parts: 



J = j ^(x) X ■ V/iQ dx = — 
Jn Jn 



J± 

q{ho) 



dx. 



(11) 



The integral J^ is interpreted as the friction factor. The 
divergence of this integral at A = is the reason for 
introducing the slip length in Eq. ^. Since, however, 
this length is very small, being measured on molecular 
scale, the integral can be evaluated by separating it into 
two parts. Near the contact line, i.e. in a ring R < r < I 
where X ^ I <^ R, the integration can be carried out 
using the linearized profile h = 9{R — r). This yields, 
asymptotically at Z ^ A, 



Ji^eirRl [9oiR-r) + 3X]-^dr: 
'r-i 



T'"IJ- <-) 



In the bulk region r < R — l, X can be neglected, and the 
integration yields, asymptotically at I <^ R, 



III. CHEMICAL SELF-PROPULSION 
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When both integrals add up, the auxiliary length I falls 
out, resulting in an expression containing the logarithm 
of the ratio of the macroscopic and microscopic scales: 



^ GttR OoR 



(14) 



An additional contour integral dependent on the un- 
known first-order function hi is contributed by the op- 
erator C when it is multiplied by ip and integrated by 
parts: 

It = - f ^{s) fc(/io) n ■ VV^/iids + (p {n ■ x) haV^hids 

- (j) (x ■ V/io) n ■ V/ii ds + (j) (n-x) V^Hq hids. (15) 

The first integral vanishes, since at ft. ^ q{h) ex h^, 
while p{h) is only logarithmically divergent. The second 
integral vanishes at /iq = as well, and so does the last 
integral, since /iq is harmonic. The remaining integral 
expresses the driving force due to the variable part of 
the contact angle 9 = 9 — Oq = —S~^n ■ Vft-i. Using 
the identity x ■ Vfto = —Oocosip, where ip is the polar 
angle counted from the direction of motion, the force is 
evaluated as 



a£a_ 
S 



J' 



(16) 



Thus, the solvability condition defining the droplet ve- 
locity reads 



J 
T 



67r 6A 



cos(?!)6i((?!))d0. (17) 



A. Surfactant distribution 



Variation of the contact angle is caused by substrate 
modification, e.g. dissolution of the surfactant adsorbed 
on the substrate in experiments of Sumino et al. |5|, Ig 
(Fig. ^ . We write the the adsorption-diffusion equation 
for the surfactant coverage on the substrate in the di- 
mensionless form 



Ct 



y\ 



H{x). 



(18) 



Here H{x) equals to 1 outside and inside the droplet 
footprint. The surfactant coverage c is scaled by the 
coverage in equilibrium with the surfactant concentra- 
tion in the continuous phase, time by the inverse adsorp- 
tion/desorption rate constant k (which are assumed for 
simplicity to be equal), and length, by yj D /k^ where D 
is the surfactant diffusivity on the substrate. 

The equation transformed to the frame moving with a 
dimensionless velocity v = U /\/ Dk along the x axis is 



VCx 



^^c-c^H{x) =0. 



(19) 



The solution of this equation can be expressed with the 
help of an appropriate Green's function through an in- 
tegral over the droplet footprint area and subsequently 
transformed into a contour integral with the help of the 
Gauss theorem llOl . For a circular contour with a dimen- 



sionless radius (Thiele modulus) a — Ryjk/D, the con- 
centration on its boundary, which determines the propa- 
gation speed, is computed in this way as 



J 



=(0) = 1 - 



27r 



— ■^va{cos 0— cos^) 



1 H sin 

4 2 



2COS^Xo|2ayi + ^sin^y^ 



Xi|2a\/l + ^sin^^ 



d?, 



(20) 



r 



where Kn (and /„ below) are modified Bessel functions. 



Assuming the contact angle to be a linear function of 



combining Eqs. ifTTjl and (jSOJ: 
M 



■ V cos 



(21) 



the surfactant coverage, 9 = —/3c, the propagation ve- where c = 1 — c and all relevant parameters, except the 
locity is computed by solving the equation obtained by droplet radius a in Eq. (|20|l . are lumped into a single 




FIG. 1: Schematic diagram of droplet motion due to modifi- 
cation of surface wettability 5^ Sj. The adsorbed surfactant 
is represented as a hydrophobic bar with a hydrophilic head. 



dimensionless combination, which can be called the mo- 
bility parameter: 



M = 



%ap 



\n 



6A ■ 



(22) 



The parameter M retains a weak logarithmic dependence 
on the droplet radius. 



B. Traveling Bifurcation 

Simplified expressions can be obtained in the limiting 
cases u <C 1 and ■;; ^ 1. The parameter i5 can be still 
small also in latter case, since the ratio of the characteris- 
tic "chemical" velocity \/Dk to the characteristic velocity 
a/r] which determines the influence of viscous stresses on 
the droplet shape is typically very small. 

The limit v -^ 1 (i.e. fast diffusion) is analogous to the 
fast inhibitor limit in reaction-diffusion systems, which is 
conducive to formation of stationary patterns |2ll | . Equa- 
tion 12U|I is expanded in this limit to the first order in v 
as 



c(0) = 3 



sm — K^ 2a sm — 
2 V 2 



1 f C 

- cos CKol2a sm - 



r 



sCr. ^_..„C 



a sin — Ki {2a sin — 
2 \ 2 



+ 0(i;2) dC, (23) 



where ( = £, — 4>. In the leading 0(1) order, the surfac- 
tant distribution is circularly symmetric. The first-order 
dipole term in Eq. H23|l is the only one contributing to the 
integral in Eq. H21|l (another term vanishing upon inte- 
gration is omitted). The angular integrals are evaluated 
using the identities 



$fc(a) 



*fe(a) 



sm 



2 k 



Kq I 2a sin ■ 



J_ d{a^k-i) 
2a 



,2fe+l 



Ki 2a sin 77 ) d0 = — 



da 
ld£fe 
2 do ' 



starting from $o(a) = '^Io{0')Ko{a). Plugging the result- 
ing expressions in Eq. (|21|l yields the condition for the 
onset of motion 



Mo' 



[h{a)K2{a)~Io{a)Ki{a)]. 



(24) 



The critical value Mq^ as a function of a is plotted in 
Fig. 121 Since the radial dependence in Eq. H24(l saturates 
when the droplet radius far exceeds the diffusional range, 
so that Mq —^ 4 at a —f 00, no droplet can move below 
this limiting value. 

Because of the logarithmic dependence of M on the 
droplet radius, the size dependence of the traveling 
threshold in Eq. (|24|) still remains implicit. An explicit 
dependence can be extracted after rewriting Eq. (I25II as 



M^fialn 



A^d 



e^ap 




In the limit a ::^ 1 when the critical value Mq in Eq. 1)24(1 
approaches the limit Mq = 4, the droplet is mobile at 
a < Xd exp(/ic;/4). This suggests that the radius of mobile 
droplets is bounded both from below and above, as it 
indeed follows from the existence of a maximum in the 
radial dependence of the critical value of the parameter 
fid, such as seen in Fig.|31 

Beyond the critical point, the velocity can be obtained 
using further terms in the expansion (|23() . The dipole 
component of the term quadratic in v vanishes, and there- 
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FIG. 2: The critical value MJ~ as a function of the dimen- 
sionless droplet radius a. Inset: blow-up near the origin on 
logarithmic scale. The droplets are mobile below this curve. 
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FIG. 3: The critical value of /i^ ^ as a function of the dimen- 
sionless droplet radius a for Ad = 10~*. 



fore there is no contribution to motion in this order. The 
dipole component of the third-order term is a^v^ where 

"3 = -^ [h{a)Ki{a) - h{a)Ko{a) - a-^h{a)K^{a)\ . 
16 

.(26) 
This coefficient is negative; hence, the bifurcation is su- 
percritical and propagation is possible at M > Mp. For 
small deviation s M2 — M — Mq > 0, the velocity is 



C. Non-diffusive Limit 

In the opposite limit when diffusion is neghgible, the 
surfactant concentration can be obtained directly by inte- 
grating along the direction of motion with the initial con- 
dition c = 1 on the advancing contact line (|(/)| < 7r/2). 
The resulting concentration on the receding contact line, 



c(0)=e"2^™"'^, 
depends only on the ratio 

a kR p 

'^ ~'v~ IT ^'5 

where p = Rkrj/j is the rescaled dimensionless droplet 
radius. The equation for r following from Eq. (|17|l is 



(27) 



(28) 



P 



r/2 



cos 6 ( 1 



-2r cos r 



T^ JO 

containing a single parameter 



r 



(29) 



M «2 



6T]Rk ^ 9qR 

m 

6A 



e'ocrp 



where 



M/i 



6 '' 



X 



P , P 
— In—, 

Ph Ah 
6A7 



kri0o 



We rewrite Eq. (|^ as 



P , P 
— In — 

Ph Ah 



tF{t). 



The function F{t) is evaluated as 
1 1 



F{t) 



[/i(2r)-L_i(2T) 



(31) 



where L„(a;) is a Struve function. The function F{t) 
increases monotonically from at r = to I/tt at r — > 00. 
The dependence of velocity on the droplet radius can 
be obtained analytically in two limiting cases correspond- 
ing to the unsaturated and saturated regimes, respec- 
tively, at small and large r. In the former case, one can 
use the approximation F{t) = r/2 + 0{t'^) to obtain 



PPh 



2ln{p/Xh) 
In the opposite limit r ^ 1, F{t) 

1 

'5 



1/2 



I/tt and 



TT p 

-In—. 
Ah 



Ph 



(32) 



(33) 



Thus, the velocity increases with droplet size in the un- 
saturated and decreases in the saturated regime, in agree- 
ment with experiment |3] and earlier computations |9|. 
For intermediate values of t, the dependence of veloc- 
ity on radius obtained by solving Eq. (j3()|l numerically is 
plotted in Fig. ^ 
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FIG. 4: The dependence of the capillary number 5 on the 
dimensionless radius p for Ah — 10"* and /ih = 1 (a) and /ih ~ 
0.1 (b). The dashed line shows the asymptotic dependence 



D. General case 



In a general case, Eq. H21() can be rewritten as M^^ — 
G{a, v) and solved after computing numerically the dou- 
ble integral 



G{a,v) 



a 



— ^va[cos(l)—cos(<t>+C,)] ^ 



^ cos(<^ + Oi^o [20^1 + ^ Sin ^ 



1 + ^smM;,, 2a^l + ^sinKi 



dC, (34) 



The function G{a, v) is plotted against v at several values 
of a in Fig. |S1 The curves peak at the ordinate at the 
bifurcation value M^^ given by Eq. (|^ . 
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additive. This leads to the equation of motion for droplet 
centers X^ 



dt 






X, - X 



k a 



IX,, - Xl 



^K,{\X,-Xk\). (36) 



This is a gradient dynamical system 



dX, 
dt 



Mj dV 
a, dXi 



evolving to minimize the potential 
1 






Xj-Xk\). 



(37) 



(38) 



In a confined region, the potential is expected to be 
minimized by a regular hexagonal pattern with spac- 
ing dependent on the number density of droplets. This 
is demonstrated by an example of evolution shown in 
Fig. |H| The confinement is effected in this computation 
by a centripetal external potential. One can see that evo- 
lution starting from a random arrangement of droplets 
evolves to regular pattern where circles mark final posi- 
tions falling on a hexagonal grid. This might be a prac- 
tical way to arrange a regular dewetting pattern on a 
homogeneous substrate. 



B. Scattering 

The surfactant depletion in the far field of a steadily 
moving droplet can be obtained by solving Eq. 119|l with c 



FIG. 5: Plots of G(a, v) defined by Eq. l|!^ as a function of 
v; the values of a axe marked at the respective curves. 



IV. DROPLET INTERACTIONS 

A. Relaxation to a Stationary Pattern 

At high diffusivities when droplets do not travel spon- 
taneously, they still can move under the influence of mu- 
tual interactions. The surfactant depletion in the far 
field of a stationary droplet (at distances far exceeding its 
radius) is well approximated by the stationary solution 
Eq. ((TH|) with c replaced by c and H{x), by the delta- 
function multiplied by the droplet area: 



c^^a'Koir). 



(35) 



The circular symmetry of the surfactant distribution 
around a single droplet is perturbed by the far field of its 
neighbors. The resulting repelling interaction induces, 
according to Eq. (|17|) . motion with the velocity propor- 
tional to the concentration gradient at the droplet lo- 
cation. If there is a number of droplets, their action is 
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FIG. 6: Trajectories of droplets moving according to Eq. I36II 
with added centripetal force. The dots mark the original ran- 
dom arrangement, and the circles, final positions. 






FIG. 7: Surfactant depletion in the far field of a stationary 
droplet (left) and droplets propagating with the speed v — 1 
(center) and v = 2 (right). The area with surfactant depletion 
above the same level is shaded, showing the depletion in the 
tail region increasing at higher speed. 
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FIG. 8: Scattering trajectories. The vertical line marks the 
location of the closest approach. 



replaced by c and H{x), by the delta- function multiplied 
by the droplet area. The solution is expressed in polar 
coordinates r, (jj centered on the droplet as 



a_^_i„^,o,^^^ / 



V 



(39) 
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FIG. 10: Dependence of the change of the scattering angle 
on velocity for F = (a) and its dependence on the separa- 
tion interval for v = 1 (b). Both curves are computed using 
Eq. (BU. 



The depletion field is strongly asymmetric, with a slower 
decay behind the droplet (Fig. [Tj) . 

Droplets moving one toward another are scattered by 
mutually repelling interaction created by the gradient of 




FIG. 9: The change of the scattering angle with distance for 
V = 1, Y = 2. The dashed line shows the scattering angle at 
infinity computed with the help of Eq. 14211 . 



the far field. The problem remains tractable in the qua- 
sistationary approximation as long as velocity induced 
by interaction is much smaller than the speed of self- 
propelled motion. Otherwise, the far field becomes de- 
pendent on the entire history of motion, and solving the 
full non-stationary problem H18I) is necessary. 

Consider as an example two droplets of equal size prop- 
agating along the x axis with identical speed v on an- 
tiparallel trajectories shifted by the interval 2Y, as in 
Fig. IHl Since a deflected droplet keeps moving on a per- 
turbed course, the scattering action is equivalent to ac- 
celeration in the direction normal to self-propelled mo- 
tion. Restricting to the quasistationary approximation, 
the dynamic equation for the deviation y normal to the 
original trajectory is therefore 



di2 



-Ma 



ffc 
dy ' 



(40) 



where the derivative of the surfactant depletion given by 
Eq. (|39() is computed at a current distance between the 
droplets at the moment t equal to 



r(t)=2 (vtf + {Y + y)' 



1/2 



(41) 



where the moment of closest approach is taken as t = 0. 
Neglecting the change of the velocity component along 
the X-axis, time t in Eq. H4()|l can be replaced by x/v. 
A typical trajectory obtained by integration is shown in 
Fig. |S1 and a more quantitative example of the change 
of the scattering angle with distance is shown in Fig. 
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Take note that, due to a faster decay of depletion ahead of 
the droplet, scattering largely accumulates already after 
the droplets have passed the point of closest approach. 

For moderate deviations, a reasonable approximation 
for the scattering angle at infinity lu — y'(oo) gives the 
formula neglecting y compared to Y: 



MYa-^ 

2v 



T 



Va;2 + y2 



Ki (7V(4 + w2)(a;2 + y2)') dx. 



(42) 



r 



Scattering angle computed with the help of this formula 
only weakly depends on velocity (Fig. llOb.'l. A much 
stronger dependence on the separation interval is shown 
in logarithmic coordinates in Fig. llUb . 



V. SURFACE FREEZING AND MELTING 

A different mechanism of spontaneous motion has been 
observed in processes of surface freezing and melting y] . 
In these experiments, liquid alkane wets silicon substrate 
better than a frozen smectic layer. Respectively, the equi- 
librium contact angle increases with growing thickness of 
smectic, and a droplet tends to slip to a lower level when 
placed at a terrace edge. 

For a droplet sitting on a terrace as shown in Fig. llir al. 
9 has distinct constant values on segments |0| < "0 and 
-0 < \4>\ < TT. The contact angle is smaller at the lower 
terrace, which we place on the right; we take this value 
as 9o and denote the contact angle at the higher terrace 
as 6*0 -I- 0. Then the velocity defined by Eq. (|17|l is 



U = Uq sin0. 



f/o 



rj 67r 



In" 



6A 



(43) 




During freezing, the change of the angular front position 
"0 due to propagation of the freezing front beneath the 
droplet with a constant velocity C obeys i?sin-0dV' — 
—C At. Due to the droplet motion, the net front velocity 
relative to the droplet center is C~U. Thus, the dynamic 
equation of the angular position is 



K-r- - Uo : — -. 

di sm ip 



(44) 



The stationary position is "00 = arcsin(C/t/o); Uq is the 
maximum speed allowing for ths equilibrium. The con- 
figuration is stable at cos0 < 0, i.e. > tt/2. Under this 
condition, the length of the line where the contact angles 
on the advancing and receding sides decreases when the 
droplet slides ahead, so that the freezing front catches 
up and the stationary configuration is restored. An- 
other equilibrium position at -0 < tt/2 is unstable. This 
is in agreement with the disposition seen in Fig. Illf b') 
taken from the experimental freezing sequence ^22j. The 
droplet travels forward, while its tail is hooked to the 
frozen terrace left behind. Stabilization at an obtuse an- 
gle, which is possible only in 3D, makes unnecessary a 
hypothetic synchronization mechanism through heat ex- 
changed invoked in earlier 2D computations |23l | . 

One could expect such an equilibrium configuration 
to be impossible during melting when the directions of 
motion of the phase transition front and the droplet slip 
given by Eq. (|17|l are opposite. This equation, however, 
is not applicable during melting transition, since, unlike 
freezing when the back part of the droplet sits on the 
terrace it has created, the droplet formed as a result of 
melting is attached to the high-energy side surface of the 
smectic layer on the melting edge, and is carried along as 
this edge propagates Q- 
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